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ABSTRACT 

We present the results of Giant Metrewave Radio Telescope (GMRT) observations to 
detect Hi in absorption towards the cores of a sample of radio galaxies. From obser- 
vations of a sample of 16 sources, we detect Hi in absorption towards the core of only 
one source, the FRII radio galaxy 3C452 which has been reported earlier by Gupta 
& Saikia (2006a) . In this paper we present the results for the remaining sources which 
have been observed to a similar optical depth as for a comparison sample of compact 
steep-spectrum (CSS) and giga-hertz peaked spectrum (GPS) sources. We also com- 
pile available information on Hi absorption towards the cores of extended radio sources 
observed with angular resolutions of a few arcsec or better. The fraction of extended 
sources with detection of Hi absorption towards their cores is significantly smaller 
(7/47) than the fraction of Hi detection towards CSS and GPS objects (28/49). For 
the cores of extended sources, there is no evidence of a significant correlation between 
Hi column density towards the cores and the largest linear size of the sources. The 
distribution of the relative velocity of the principal absorbing component towards the 
cores of extended sources is not significantly different from that of the CSS and GPS 
objects. However, a few of the CSS and GPS objects have blue-shifted components 
>1000 km s _1 , possibly due to jet-cloud interactions. With the small number of de- 
tections towards cores, the difference in detection rate between FRI (4/32) and FRII 
(3/15) sources is within the statistical uncertainties. 

Key words: galaxies: active - galaxies: evolution - galaxies: nuclei - galaxies: jets - 
radio lines: galaxies 
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1 INTRODUCTION 

It is widely believed that the source of energy of active galac- 
tic nuclei (AGN) is the accretion of matter onto a super- 
massive black hole with an accretion disk, located in the 
centre of the galaxy. While there have been many theoret- 
ical and numerical studies of accretion flows (e.g. Hawley 
2011), observational studies of the gaseous components in 
the circumnuclear regions could provide useful insights to- 
wards understanding the dynamics and properties of the gas 
that might be fuelling the AGN. At radio frequencies an im- 
portant method of probing this region is via 21-cm Hi ab- 
sorption towards compact radio sources, such as the cores of 
radio galaxies and quasars, and the compact steep-spectrum 
(CSS) and gigahertz peaked-spectrum (GPS) sources. The 
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CSS and GPS objects (O' Dea 1998) typically have sizes 
of <10 and 1 kpc respectively, and are believed to be the 
young (<10 yr) progenitors of the larger radio galaxies and 
quasars which extend up to a few Mpc, with the largest 
sources being typically ^10 8 yr old (Jamrozy et al. 2008; 
Konar et al. 2008). Such studies could also help us study 
the evolution of the gaseous component with the age of the 
radio source and also test consistency of these properties 
with the unified schemes for AGN (Pihlstrom, Conway & 
Vermeulen 2003; Gupta & Saikia 2006b; Curran & Whiting 
2010). 

Studies of Hi absorption towards radio galaxies and 
quasars have shown that the absorption lines are more of- 
ten detected towards CSS and GPS objects, with previous 
studies (Pihlstrom, Conway & Vermeulen 2003; Gupta et 
al. 2006; Chandola, Sirothia & Saikia 2011) showing an ap- 
parent anti-correlation between the integrated optical depth 
and source size. For an assumed fixed value of the covering 
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Table 1. Radio sources observed with the GMRT. 



Source 


Alt. 


Redshift 


Si. 4 


Pi. 4 


Type 


LAS 


LLS 


LAS 


VLBI 


name 


name 




mjy 
(mjy) 


xlO 25 
(WHz -1 ) 




(") 


(kpc) 


ref. 


ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J0209+3547 


4C+35.03 


0.0377 


2087 


0.70 


FRI 


83 


61 





1 


J0223+4259 


3C 66B 


0.0213 


6217 


0.65 


FRI 


678 


288 


1 


2 


J0313+4120 


S4 


0.1340 


487 


2.38 


FRII 


570 


1340 


2 


3 


J0334+3921 


4C+39.12 


0.0206 


1183 


0.12 


FRI 


402 


165 


1 


2 


J0418+3801 


3C111 


0.0485 


15121 


8.59 


FRII 


195 


182 





4 


J0656+4237 


4C42.22 


0.0590 


937 


0.80 


FRI 


67 


76 







J0748+5548 


DA240 


0.0357 


2237 


0.67 


FRII 


2040 


1428 


3 


5 


J0758+3747 


3C189 


0.0428 


2637 


1.16 


FRI 


144 


120 





2 


J1147+3501 




0.0631 


815 


0.80 


FRI 


780 


936 


4 


2 


J1628+3933 


3C338 


0.0304 


3595 


0.78 


FRI 


112 


67 





2 


J1744+5542 


NGC6454 


0.0312 


810 


0.18 


FRI 


200 


123 


5 


3 


J1835+3241 


3C 382 


0.0579 


5548 


4.55 


FRII 


171 


189 


6 


2 


J1836+1942 


PKS 


0.0161 


1458 


0.09 


FRI 


316 


102 


1 




J1842+7946 


3C390.3 


0.0555 


11537 


8.73 


FRII 


235 


250 





6 


J2245+3941 


3C 452 


0.0811 


10383 


17.25 


FRII 


256 


386 


7 


2 


J2338+2701 


3C465 


0.0302 


7590 


1.63 


FRI 


487 


290 


1 


7 



Column 1: source name; column 2: alternative name; column 3: redshift ; column 4: flux density at 1.4 GHz in mjy; column 5: 

luminosity at 1.4 GHz; column 6: radio morphology i.e. FRI or FRII; columns 7 and 8: largest projected angular (LAS) and linear 

(LLS) size in arcsec and kpc respectively, where LAS has been estimated from the outer hotspots in the FRII sources, and from the 

outermost contours in the FRI sources; column 9: references for the LAS; column 10: references for the VLBI structure. 

References for the largest angular sizes - 0: This paper; 1: NRAO VLA Sky Survey (NVSS; Condon et al. 1998); 2: de Bruyn (1989); 3: 

Willis, Strom & Wilson (1974); 4: Giovannini et al. (1999); 5: Bridle & Fomalont (1978); 6: Leahy & Perley (1991); 7: Gupta & Saikia 

(2006a) 

All red-shift information is from NED (NASA Extragalactic Database) except for J1744+552 and J1842+7946. The redshift for 
J1744+552 has been taken from Bridle & Fomalont (1978) and for J1842+7946 has been taken from Eracleous & Halpern (2004). 
Flux densities at 1.4 GHz are determined using the maps from NVSS ; Luminosities have been calculated from these flux densities 
assuming a spectral index of 1. 

References for the VLBI structure: 1: Lara et al. (1997); 2: Giovannini et al. (2001); 3: Britzen et al. (2007); 4: Linfield (1987); 5: 
Saripalli et al. (1997); 6: Alef et al. (1996); 7: Venturi et al. (1995) 

et al. (2001), and more recently Emonts et al. (2010). To 
understand any evolution in the Hi component of the cir- 
cumnuclear gas as the radio source grows from the GPS and 
CSS scales to larger dimensions, it is necessary to probe 
the gas on similar scales from radio observations probing 
the compact regions with comparable optical depths. This 
can be achieved via 21-cm Hi observations towards the ra- 
dio cores of large radio sources with resolutions of about a 
few arcsec or better so that the core is resolved from the 
more extended emission and can be unambiguously identi- 
fied. Such observations will also enable us to determine any 
differences in the Hi properties of FRI and FRII sources. 
There have been suggestions that the mode of fuelling the 
nucleus and the structure of the torus/disk may be different 
for the FRI and FRII radio sources, with those in FRIs 
being geometrically thin (e.g. Morganti et al. 2001). 

In this paper we present the results of 21-cm Hi ab- 
sorption observations towards the cores of FRI and FRII 
radio galaxies with the Giant Metrewave Radio Telescope 
(GMRT), and combine our results with those in the liter- 
ature to examine any evolution in the Hi properties in the 
circumnuclear region with source size. Possible differences 
between FRI and FRII radio galaxies, and the distribu- 
tion of the velocities of the absorbing components relative 
to the systemic velocities have also been examined. The ob- 
servations are described in Section [2] The discovery of 21- 



fraction and spin temperature this corresponds to an anti- 
correlation with the neutral hydrogen column density. Gupta 
& Saikia (2006b) found the 21-cm absorption detection rate 
to be higher towards the CSS and GPS sources identified 
with galaxies than with quasars. They find that there is a 
tendency for the detection rate as well as the column den- 
sity for galaxies to increase with core prominence, a statis- 
tical indicator of the orientation of the jet axis to the line 
of sight. This can be understood in a scenario where radio 
sources are larger than the scale of the circumnuclear Hi disk 
so that the lines of sight to the lobes at very large inclina- 
tions do not intersect the disc. On the other hand, Curran & 
Whiting (2010) and more recently Allison et al. (2012) have 
investigated the effect of UV luminosity on the 21-cm ab- 
sorption detection rate. They found that after excluding the 
sources with UV luminosities, Luv >10 23 WHz" 1 , the 21- 
cm absorption detection rate for compact radio sources is 
only marginally higher than for other radio sources. Curran 
& Whiting suggest that higher 21-cm absorption detection 
rate amongst compact objects is probably due to their gen- 
erally low UV luminosity rather than their compact sizes or 
orientation. 

Compared with the studies of Hi absorption towards 
CSS and GPS objects, the extended radio sources have re- 
ceived less attention. However, there have been a few no- 
table exceptions like van Gorkom et al. (1989), Morganti 
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cm absorption towards the core of the FRII radio galaxy 
J2245+3941 (3C452), which is a part of our sample, has 
been reported earlier (Gupta & Saikia 2006a). In this pa- 
per we present the results for the remaining fifteen sources. 
In Section [3] we combine our results with those of similar 
searches and discuss the statistical trends in gas properties 
with respect to radio-source characteristics. The results are 
summarized in Section [4] 



2 OBSERVATIONS 

Table [1] lists the 16 FRI and FRII radio galaxies which we 
observed with the GMRT over the period of 2004-2010. 
Table Q] also lists some of their observed properties. With 
the available observing time, these 16 were selected largely 
from the 3CR and B2 samples to have a core flux den- 
sity at 1400 MHz >100 mjy with resolutions of a few arc- 
sec or better so that the core could be distinguished from 
the more extended bridge/jet emission. The 16 sources ob- 
served with the GMRT consist of 10 FRI and 6 FRII radio 
sources. Their projected linear sizes range from ~ 61 kpc to 
as large as ~ 1.4 Mpc (Ho =71 km s" 1 Mpc~\ Q m =0.27, 
n vac =0.73; Spergel et al. 2003). 

The observations were made with a baseband band- 
width (BB BW) of 4 MHz (velocity range ~ 800 km s" 1 ) 
split into 128 channels (resolution ~ 7 km s _1 ) , except for 
J0313+4120 for which a BB BW of 8 MHz (resolution ~ 

15 km s" 1 ) was used (see Table[2]). For the observations of 
J2245+3941 (3C452) on 2005 Dec 11 the then new high- 
resolution mode (resolution ~ 3.6 km s _1 ) of the GMRT 
correlator was used (cf. Gupta & Saikia 2006a). The ex- 
pected redshifted 21-cm frequency was placed at the centre 
of the BB BW by tuning the local oscillator chain. The stan- 
dard flux density and bandpass calibrators (3C 48, 3C 147 
and 3C 286) were usually observed every three hours to cor- 
rect for variations in amplitude and bandpass. Observations 
of the target source were preceded and followed by observa- 
tions of a phase calibrator, a compact radio source, typically 
every 35 min. 

The Astronomical Image Processing System (MPS) 
package was used to reduce the data. After flagging bad 
data, including those affected by Radio Frequency Interfer- 
ence (RFI), and calibrating the data, a continuum image 
of the source was made by averaging ~60 line-free chan- 
nels. The images were made with the highest resolution, 
which was typically ~2 arcsec with the GMRT at L-band, 
because of our interest in detecting absorption towards the 
compact components such as the radio cores. The data were 
self-calibrated to produce a satisfactory continuum image; 
after which the complex gains were applied to all the fre- 
quency channels and the continuum image was subtracted 
from the visibility data cube. Spectra were extracted at the 
pixels of maximum brightness in the cores and the hot-spots. 
This was done independently for both Stokes RR and LL 
to check for consistency. The two polarization channels were 
later combined to produce the final Stokes I spectrum which 
was shifted to the heliocentric frame. The results for all the 

16 sources, including the radio galaxy J2245+3941 (3C 452) 
are summarised in Tabled 

The final radio continuum images, i.e. the ones used for 
continuum subtraction, and the 21-cm Hi absorption spectra 
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Figure 1. Hi column density distribution of the 'cores sample'. 
Shaded areas represent the detections while the unshaded areas 
represent 3<r upper limits. 
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Figure 2. Upper panel: Hi column density distribution of the 
'CSS & GPS sources' from Gupta et al. (2006). Lower panel: Hi 
column density distribution of CSS & GPS sources with N(H i) 
or upper limits to it greater than 1.5xl0 20 cm" 2 . Shaded areas 
represent the detections while the unshaded areas represent 3it 
upper limits. 
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Table 2. Observational details and results of the search for associated Hi absorption. 



Source Date Obs. freq BW t; ut Beam Comp. Peak flux cr r AT(Hl) 

(MHz) (MHz) (hrs) ("x" °) (mjy/b) (mjy/b) (lO^ 3 ) (10 20 c m - 2 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



J0209+3547 


2006 


Jan 


1368.76 


4 


3.5 


2.41x1.91 


+20 


Core 


102 





.91 


<26.7 


<5.15 


J0223+4259 


2006 


Jan 


1390.84 


4 


1.6 


2.64x2.25 


—83 


Core 


170 


1 


.10 


<19.5 


<3.77 


J0313+4120 


2004 


Jan 


1252.56 


8 


3.0 


2.62x2.15 


—53 


Core 


383 


1 


.12 


<8.7 


<1.68 


J0334+3921 


2006 


Jan 


1391.75 


4 


4.9 


2.70x1.99 


+53 


Core 


304 





.85 


<8.4 


<1.62 


J0418+3801 


2006 


Jan 


1354.70 


1 


5.2 


3.15x2.38 


+34 


EHS 


1393 


0. 


.90 


<2.1 


<0.41 


















Core 


1103 


1. 


.05 


<3.0 


<0.58 


















WHS 


250 


1 


.03 


<12.3 


<2.37 


J0656+4237 


2010 


Jan 


1341.27 


4 


4.9 


3.63x2.80 


+64 


Core 


224 


0. 


.53 


<6.9 


<2.95 


J0748+5548 


2010 


Jan 


1371.49 


4 


4.4 


2.70x2.16 


+25 


Core 


209 


0. 


.38 


<5.4 


<2.78 


J0758+3747 


2006 


Jan 


1362.06 


4 


5.2 


3.12x2.07 


+43 


Core 


184 


1. 


.01 


<16.5 


<3.19 


J1147+3501 


2010 


Jan 


1336.06 


4 


4.4 


3.69x2.38 


+38 


Core 


351 





.57 


<4.8 


<1.96 


J1628+3933 


2010 


Jan 


1378.56 


4 


4.8 


3.42x2.17 


+61 


EP 


67.6 





.93 


<41.4 


<7.99 


















Core 


179 


1. 


.14 


<19.2 


<3.70 


















WP 


69.5 


0. 


.76 


<32.7 


<6.32 


J1744+5542 


2010 


Jan 


1377.43 


4 


4.0 


4.24x3.18 


-75 


Core 


313 


1 


.18 


<11.4 


<2.20 


J1835+3241 


2006 


Mar 


1342.70 


4 


4.5 


2.68x2.10 


+35 


EHS 


203 


1 


.05 


<15.6 


<3.01 


















Core 


240 


1. 


.17 


<14.7 


<2.84 


















WHS 


37 


1 


.12 


<90 


<17.4 


J1836+1942 


2004 


Jan 


1397.84 


4 


2.0 


3.04x2.10 


-70 


Core 


224 


2. 


.10 


<28.2 


<5.44 


J1842+7946 


2010 


Jan 


1345.72 


1 


4.9 


3.88x2.32 


+17 


NHS 


177 


1 


.02 


<17.4 


<3.34 


















Core 


165 


1. 


.15 


<18.9 


<4.02 


















SHS 


839 


1. 


.90 


<6.9 


<1.31 


J2245+3941* 


2005 


Dec 


1313.85 


4 


5.6 


3.14x2.25 


+52 


EHS 


94 


1. 


.10 


<35.1 


<6.77 


















Core 


194 


1 


.00 


58.0 


6.39 ±0.13 


















WHS 


128 


1 


.10 


<25.8 


<4.98 


J2338+2701 


2004 


Jan 


1378.74 


4 


1.0 


2.89x2.24 


-50 


Core 


302 


1. 


.63 


<16.2 


<3.13 



Column 1: source name; column 2: year and month of observations; column 3: redshifted 21-cm frequency; column 4: baseband 
bandwidth (BW); column 5: observing time in hr (excluding calibration overheads); columns 6, 7 and 8: restoring beam, radio 
component (EHS = eastern hotspot, WHS = western hotspot, NHS = northern hotspot, SHS = southern hotspot, EP = eastern peak, 
WP = western peak and Core= radio core) and peak brightness in mjy/beam for the continuum image made using line-free channels 
respectively; column 9: rms noise in the spectrum; column 10: optical depth estimate, and column 11: Hi column density in units of 
10 20 cm -2 , assuming T s =100 K and a covering factor, f c of unity; upper limits are 3a values, assuming Av=100 km s — 1 . 
* The 21-cm absorption towards J2245+3941 (3C452) has velocity of —2.8 km s" 1 relative to its optical systemic velocity. 



for the 15 galaxies with non-detections are presented in the 
Appendix (Figs. lA"lT - |A15[) . The large-scale structure is seen 
in most of the images except for J0313+4120, J0748+5548, 
J1147+3501, J1744+5542 and J1836+1942 where the diffuse 
emission is largely resolved out in our GMRT images. For 
all sources the rms noise in the spectra, the column density 
for the lone detection, 3C452, and 3-<r upper limits to the 
optical depth towards all the other cores, and also towards 
the hotspots or peaks in 5 of the 16 sources are presented in 
Table. [2] The Hi column density, iV(fil), has been estimated 
using the relation 

,„ T s f r(v)dv 
iV(Hi) = 1.823 x 10 1S J , cm" 2 (1) 

Jc 

and assuming T s =100 K and / c =1.0, where T B , r and f c 
are the spin temperature, optical depth at a velocity v and 
the fraction of background emission covered by the absorber 
respectively. In this paper wherever we mention Hi column 
density, it refers to the pseudo Hi column density with the 
above assumptions. The upper limits to iV(Hl) are 3-<r values 



and have been estimated, as stated above, assuming T a =100 
K, Au=100 km s -1 and / o =1.0. 



3 DISCUSSION 

In order to estimate reliably the fraction of sources with 
21-cm Hi detection, it is important to consider surveys of 
sources which report the detections as well as the non- 
detections. In our sample of 16 radio galaxies observed with 
the GMRT we detect 21-cm absorption associated with the 
core of only 1 radio galaxy i.e. J2245+3941 (3C 452) (Gupta 
& Saikia 2006a). To improve upon the statistics, we have 
considered the observations of large (>15 kpc) sources by 
van Gorkom et al. (1989), Morganti et al. (2001) and Emonts 
et al. (2010) who have reported both the detections and 
non-detections in their samples, and have summarized their 
results in Table [3] which has a total of 31 sources. We have 
included only those sources which have been observed with 
few arcsec resolution (corresponding to a beam size of <15 
kpc for our red-shift range) in order to minimise the effect of 
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Table 3. Sources included in 'cores sample' from literature 



Source 


Alt. 


Redshift 


Si. 4 


Pl.4 


LAS 


LLS 


LAS 


VLBI 


Type 


AT (Hi) 


V shift 
(km s" 1 ) 


Ref. 


name 


name 






(mJy) 


xlO 25 
(WHz" 1 ) 


(") 


(kpc) 


ref. 


ref. 




xlO 20 
(cm -2 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


juuoi U-Luy 


3C15 





.0730 


4092 


5.46 


33.7 


46 






FRII 


<0.50 




M 


100^4-3091 


NOr i 31 fi 

IN V_x vjOIJ 





.0165 


2371 


0.15 


3180 


1152 


2 


\ 


FRI 


7 Q 


462 


\< 


100^7—01 93 
j \jyjo 1 u j. zo 


3C29 





.0450 


5354 


2.60 


139 


121 


I 




FRI 


<9.39 




M 


101 07-4-3994 


3C31 





.0169 


4757 


31 


900 


305 


3 




FRI 


<0 6 






101 0f\ 0190 


3P40 





.0180 


3870 
oo ( y 


90 


330 


122 






FRI 








1030&-J-040fi 


3P7S 





.0287 


73fi0 

/ ouu 


1 42 


210 


"1 1 

1 L.J 




9 


FRI 


<^0 73 




M 


10390-4-3047 







.0243 


1419 


90 




J-OU 


u 


Q 
o 


FRI 






\? 


1D331 -t-0933 


3C88 





.0302 


5030 


1.08 


191 


114 






FRII 


<5.06 




M 









.0399 


^830* 
ooou 


2 21 


97fi 


215 


u 




FRI 


^1 1 A p. 




M 


10^1 Q— 4^4fS 


PicA 





.0351 


fifiOOO* 
uuuuu 


19.24 


458 


316 


7 


4 


FRII 


<2.43 




M 


10^91 ^941 


MRC 





.0511 


3400* 


2 15 


47 7 


-] 7 







FRI 


<r^3 fiO 




M 


10097-1-90^^ 







.0253 


297 


0.04 


747 


376 






FRI 


<69 






11 00(^-4-34^3 


3P930 





.1005 


4645 


12.07 


2307 


4221 


9 


g 


FRII 


3 7 


162 


Q 


11 1 1 1 _L9fic;7 


Nnr , 3 c >fi3 

IN V_t uouuo 





.0331 


99 


0.03 


72 


47 


10 




FRI 


^"3 7 






11 1 4^4-1 Q3fi 


3C264 





.0217 


5494 


0.60 


510 


234 


11 


7 


FRI 


<1.54 




Q 


11 9^0-1-1 993 


M87 





.0044 


1 ^3490 


0.67 


896 


81 


12 


8 9 


FRI 


<1.35 




Q 


11 9^4—1 933 


3C278 





.0150 


7622 


0.39 


128 


38.7 


13 


10 11 


FRI 


<12.32 




M 


11 ^91 _4349 


NOCSOOO 





.0114 


5860* 


0.17 


770 


177 


9 


5 


FRI 


2.70 


4 


M 


J1323+3133 


NGC5127 





.0162 


741 


0.04 


688 


223.6 


10 


1 


FRI 


4.9 


33 


E 


J1324+3622 


NGC5141 





.0174 


886 


0.06 


52 


18.15 


11 




FRI 


3.0 


134 


E 


J1336-3356 


IC4296 





.0125 


11960* 


0.43 


1869 


470.9 


1 


5 


FRI 


<7.43 




M 


J1516+0701 


3C317 





.0345 


5608 


1.57 


41.4 


28.0 


1 


5 


FRI 


<1.74 




G 


J1632+8232 


NGC6251 





.0247 


2201 


0.31 


3150 


1547 


15 


12,13 


FRI 


<2.51 




G 


J1644-7715 


MRC 





.0427 


6500* 


2.79 


259.8 


215.9 


1 




FRII 


<8.63 




M 


J1720-0058 


3C353 





.0304 


54835 


11.84 


252.9 


151.7 


1 


12 


FRII 


42 


-100 


M 


J1806+6949 


3C371 





.0510 


2138 


1.35 


280 


275 


16 


14 


FRI 


<2.90 




G 


J1952+0230 


3C403 





.0590 


5855 


5.00 


98.5 


111 


17 




FRII 


<17.82 




M 


J2118+2626 


NGC7052 





.0156 


212 


0.01 


149 


46.6 


18 




FRI 


<2.4 




E 


J2157-6941 


PKS 





.0283 


30390* 


5.66 


78.8 


47.3 


19 


15 


FRII 


<1.90 




M 


J2223-0206 


3C445 





.0559 


5405 


4.10 


570.0 


611.0 


20 




FRII 


<8.06 




M 


J2231+3921 


3C449 





.0171 


2847 


0.19 


1500 


514.5 


3 




FRI 


<1.2 




E 



Column 1: source name; column 2: alternative name; column 3: redshift; column 4: 1.4GHz flux density column 5: 1.4GHz luminosity in 
the rest frame of the source; columns 6 & 7: largest projected angular (LAS) and linear (LLS) size in arcsec and kpc respectively; 
column 8: references for the LAS; column 9: references for the VLBI structure; column 10: radio morphology classification; column 11: 
Hi column density or a 3cr upper limit to it; column 12: the shift of the primary Hi component relative to the systemic velocity as 
measured from the optical emission lines, with a negative sign indicating a blue-shift, and column 13: references for the Hi observations. 
References for the LAS: EMorganti et al. (1993); 2: Bridle et al. (1976); 3: Andernach et al. (1992); 4: Saikia et al. (1986); 5: Bridle et 
al. (1991); 6: Jones & McAdam (1992); 7: Perley, Roser & Meisenheimer (1997); 8: Ekers et al. (1981); 9: Mack et al. (1997); 10: Parma 
et al. (1986); 11: Gavazzi, Perola & Jaffe (1981); 12: Owen, Eilek & Kassim (2000); 13: Baum et al. (1988); 14: Fanti et al. (1986); 15: 
Perley, Bridle & Willis (1984); 16: Wrobel & Lind (1990); 17: Dennett-Thorpe et al. (1999); 18: Morganti et al. (1987); 19: Fosbury et 
al. (1998); 20: Leahy et al. (1997) 

References for Hi observations: M: Morganti et al. (2001); G: van Gorkom et al. (1989); E: Emonts et al. (2010) 

Flux densities at 1.4 GHz are determined using maps from NVSS except the sources with *, values for which are from Wall & Peacock 
(1985); Luminosities have been calculated from these flux densities assuming spectral index of 1. 

References for the VLBI structure: 1: Giovannini et al. (2001) ; 2: Jones (1984); 3: Giovannini et al. (2005); 4: Tingay et al. (2000); 5: 
Venturi et al. (2000); 6: Schilizzi et al. (2001); 7: Lara et al. (1997); 8: Chang et al. (2010); 9: Kovalev et al. (2007); 10: Marecki et al. 
(2003); 11: Bondi et al. (2004); 12: Fomalont et al. (2000); 13: Jones et al. (1986); 14: Gabuzda et al. (2004); 15: Ojha et al. (2004) 



large scale structures on our estimates. We have also avoided 
inclusion of QSOs, BL Lacs, Spirals or Seyferts in the sam- 
ple as we wanted to have a sample of radio galaxies with 
their hosts being ellipticals. Thus combining our observa- 
tions of 16 sources listed in Table[T]and the 31 sources listed 
in Table[3] we obtain a sample of 47 radio sources, called the 
'cores sample'. We also list 9 sources in Table [4] for which Hi 
21 cm absorption has been reported towards their core from 
observations with angular resolution of a few arcsec but not 



included in our 'cores sample' for statistical analysis since 
in these cases non-detections were not published. 



3.1 Frequency of occurrence 

A comparison of the results of Hi 21-cm absorption observa- 
tions towards the cores of radio galaxies with those obtained 
for CSS and GPS objects could help explore possible evo- 
lution of the content of Hi gas in the host galaxies of radio 
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Table 4. Sources reported in literature for cores with Hi detection but not included in our 'cores sample' 



Source 


Alt. 


Red-shift 


LAS 


LLS 


LAS 


Type 


JV(Hi) 


V 'shift 


Ref. 


name 


name 




(") 


(kpc) 


Ref. 




xlO 20 
(cm -2 ) 


(kms -1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



J0319+4130 


NGC1275 


0.0176 


198 


70 


1 


FRI 


2.34 


2797 


YRS 


J0840+2949 


4C29.30 


0.0647 


520 


639 


2 


FRI 


47 


84 


CSG 


J0918-1205 


Hydra A 


0.0549 


310 


326 


3 


FRI 


42 


-180 


D 


J1219+0549 


NGC4261 


0.0075 


492 


75 


4 


FRI 


7.2 


-3 


J 


J1247+6723 




0.1073 


671 


1301 


5 


FRII 


6.73 


57 


S 


J1325-4301 


Cen A 


0.0018 


30564 


563 


6 


FRI 


50.1 


6 


VGH 


J1352+3126 


3C293 


0.0450 


216 


189 


7 


FRI 


78.9 


3 


B 


J1531+2404 


3C321 


0.0961 


286 


503 


8 


FRII 


92.3 


235 


CSDN 


J1959+4044 


Cygnus A 


0.0561 


127 


137 


9 


FRII 


25.4 


165 


C 



Column 1: source name; column 2: alternative name; column 3: redshift; columns 4 & 5: largest projected angular (LAS) and linear 
(LLS) size in arcsec and kpc respectively; column 6 : references for the LAS; column 7: Hi column density; column 8: radio morphology 
classification; column 8: the shift of the primary Hi component relative to the systemic velocity as measured from the optical emission 
lines, with a negative sign indicating a blue-shift, and column 10: references for the Hi observations. 

References for the structural information and LLS: 1: Pedlar et al. (1990); 2: Chandola, Saikia & Gupta (2010) 3: Ekers & Simkin 
(1983); 4: Condon & Brodcrick (1988); 5: Lara et al. 2001 ; 6: Junkes et al. (1993); 7: Bridle, Fomalont & Cornwell (1981); 8: Chandola 
ct al. (2012); 9: Pcrley, Drchcr & Cowan (1984). 

References for Hi observations: D: Dwarakanath, Owen & van Gorkom (1995); J: Jaffe & McNamara (1994); S: Saikia, Gupta & Konar 
(2007); C: Conway & Blanco (1995); YRS: De Young, Roberts & Saslaw (1973); VGH: van der Hulst, Golisch & Haschik (1983); 
CSDN: Chandola et al. (2012); CSG: Chandola, Saikia & Gupta (2010); B: Bcswick, Pedlar & Holloway (2002) 




Log Linear Size (kpc) 

Figure 3. Hi column density as a function of projected linear 
source size. Arrows mark the 3<r upper limit to the column den- 
sity. The unfilled symbols are FRI sources while filled symbols 
are FRII sources. Circles denote the sources from our GMRT 
observations while squares, triangles with vertex up and trian- 
gles with vertex down denote sources from Emonts ct al. (2010), 
van Gorkom et al. (1989) and Morganti et al. (2001) respectively. 
Rhombuses are sources which are not included in our cores sam- 
ple but have detection towards the cores of larger radio sources 
with ~ arcsec resolution. 



sources as the sources evolve. The column densities or upper 
limits to these in our 'cores sample' range from ~0.5xl0 20 
cm~ 2 to 69xl0 20 cm~ 2 . The median column density sen- 
sitivity of our 'cores sample', including the upper limits is 
3.1 xlO 20 cm~ 2 . For comparison, we consider the compila- 
tion of CSS and GPS sources from Gupta et al. (2006). From 
this compilation of 62 CSS and GPS sources we consider 60 
sources. We exclude J1415+1320 and J1945+7055 as these 
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Figure 4. The distribution of velocities for the principal Hi ab- 
sorption components with respect to the systemic velocity esti- 
mated from optical lines. Ones with bold boxes represent the dis- 
tribution for CSS & GPS sources, while the shaded areas denote 
the cores of FRI and FRII radio sources. In the cores, Perseus A 
(NGC1275) where the principal relative velocity 2797 km/s, has 
been put in bin at 2000 km/s. 



were reported as individual detections. To enlarge this sam- 
ple further we also consider 21-cm absorption searches based 
on samples of steep-spectrum sources reported in the liter- 
ature since 2006. We consider measurements of 17 sources 
from Chandola, Sirothia & Saikia (2011), 5 compact sources 
(B0258+35, B0648+27, B1122+39, B1447+27, B1557+26) 
with steep spectrum from Emonts et al. (2010), and 2 
sources from Carilli et al. (2007). The measurements of 10 
compact sources with steep spectrum are available from Al- 
lison et al. (2012) but detailed size information is not avail- 
able for these. So we are not able to include these in our 
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sample. Thus we have a sample of 84 CSS and GPS sources. 
This sample has column densities or upper limits to these 
ranging from ~0.06xl0 20 cm" 2 to 125 xlO 20 cm" 2 with a 
median value of 1.9xl0 20 cm~ 2 assuming that the values for 
the non-detections are close to the upper limits. To ensure 
that the sources are observed to similar sensitivity limits 
so that the distributions of column densities are similar for 
both the 'cores' and the 'CSS and GPS' samples, we con- 
sidered all the sources from 'CSS and GPS sample' with a 
column density estimate or an upper limit to it which is at 
least 1.5xl0 20 cm" 2 . This leaves us with a sample of 49 CSS 
and GPS sources with a median column density sensitivity 
of 3.8xl0 20 cm -2 , assuming again that the values for the 
non-detections are close to the upper limits. 

The detection rate for the 'cores sample' is rather low 
(7/47; ~15%; see Fig. [1]) compared with the detection rate 
for compact CSS and GPS sources (28/49; ~57%; see Fig. [2] 
lower panel). Considering the entire sample of 84 CSS and 
GPS objects, the detection rate (31/84; ~37%; see Fig.rjup- 
per panel ) is again significantly higher than for the 'cores 
sample'. All but one of the sources from 'CSS and GPS 
sample' with column density values <1.5xl0 20 cm -2 repre- 
sent sensitive observations of CSS and GPS objects with no 
detection of Hi in absorption. We have also examined the 
VLBI-scale structure (references given in Table [1] and Ta- 
ble [3]) of the cores of the larger sources ('cores sample') and 
find that they tend to have a core-jet structure where most 
of the flux density is from the VLBFscale core component 
which is usually smaller than the CSS and GPS sources sug- 
gesting that the difference in Hi absorption detection rate 
is not due to different covering factors. Although it would 
be useful to confirm it from a larger sample, clearly the de- 
tection rate of Hi in absorption towards the cores of larger 
sources is smaller than for the CSS and GPS objects, sug- 
gesting an evolution in the gaseous content of the host galax- 
ies . 

It is also interesting to examine the Hi detection rate 
for FR I and FR II radio sources in case this reflects any dif- 
ferences in either the torus/disk or fuelling processes (e.g. 
Morganti et al. 2001; Emonts et al. 2010). In the 'cores sam- 
ple' Hi is detected in absorption towards 4 out of 32 (~13%) 
FR I objects, compared with 3 out of 15 (~20%) for the FRII 
sources. Within statistical errors, there is no significant dif- 
ference in the detection rate towards the cores for FRI and 
FRII radio sources. These have been observed with similar 
sensitivity. 

3.2 Column density vs size 

As mentioned earlier, the Hi column density in CSS and 
GPS objects appears to be anti-correlated with source size 
( Pihlstrom, Conway & Vermeulen 2003; Gupta et al. 2006; 
Chandola, Sirothia & Saikia 2011), but shows no evidence 
of any dependence on either redshift or luminosity (Gupta 
et al. 2006). On sub-galactic scales, the observed anti- 
correlation has been interpreted by Pihlstrom, Conway & 
Vermeulen (2003) to be due to gas distributions with a ra- 
dial power law density profile. It is also relevant to note here 
that Emonts et al. (2010) observed a sample of radio galax- 
ies to detect Hi in emission and found an inverse relation 
between Hi mass and sources size. 

In Fig. [3] we plot the 21-cm Hi column density towards 



the cores against the projected linear size in kpc for the 
sample of 47 radio galaxies. In this figure, we also plot an 
additional nine sources where Hi detection has been reported 
from observations with angular resolutions of a few arcsec for 
single sources (see Table [3}, giving a total of 16 detections. 
The linear sizes range from tens of kpc to over a Mpc while 
the column densities vary from 0.5 x 10 20 cm -2 to 69xl0 20 
cm -2 . There does not seem to be any significant relation 
between the Hi column density towards the cores and the 
largest projected linear sizes of the sources. 

3.3 Relative velocity of absorbing gas 

The relative velocity of the absorbing gas could be signifi- 
cantly blue shifted if it is due to gas which has been accel- 
erated by interaction with the radio jet, while infalling ma- 
terial fuelling the central engine would appear red shifted. 
An important component of the unification scheme for FR II 
radio galaxies is the presence of a circumnuclear disk/torus. 
Hi gas associated with this and rotating around the nucleus 
may appear both blue and red shifted relative to the sys- 
temic velocity. It is worth noting that CSS and GPS objects 
often show complex line profiles with both blue- and red- 
shifted Hi absorption components. 

We compare the relative velocity of the principal ab- 
sorbing component for the 16 detections towards cores of 
larger sources with the 33 detections(including 2 sources 
from Gupta et al. (2006) which we left out earlier while do- 
ing statistical analysis of detection rates) towards CSS and 
GPS objects. While making these comparisons it is also rel- 
evant to note that different ionization lines may yield some- 
what different systemic velocities. For example in 3C 452, 
the low-ionization forbidden emission lines are broader and 
more red-shifted than the high-ionization ones. The typical 
error in optical systemic velocity estimate for the sample is 
~ 90 km s _1 . Fig. [4] shows that although a few CSS and 
GPS objects have large blue-shifted velocities of >1000 km 
s" 1 possibly due to jet-cloud interactions, the distributions 
for the cores and CSS and GPS objects are not significantly 
different. Also one of the objects with the highly blue-shifted 
component, J1815+6127, has a large error of ~900 km s _1 
(Vermeulen & Taylor 1995) in its optical systemic veloc- 
ity. The highly red-shifted component seen towards the core 
of NGC1275 in the Perseus cluster could be due to a gas 
cloud or galaxy in the intracluster medium moving towards 
NGC1275 (De Young, Roberts & Saslaw 1973). 



4 SUMMARY 

We have presented the results from a search for 21-cm ab- 
sorption towards the cores of 16 FR I and FR II radio galax- 
ies. From our search we have reported earlier one new detec- 
tion, i.e. towards the core of the radio galaxy J2245+3941 
(3C452) (Gupta & Saikia 2006a), and have now presented 
the upper limits on N(Hl) for the remaining 15 sources. We 
have combined our results with those of similar 21-cm ab- 
sorption searches in the literature to obtain a larger sample, 
called the 'cores sample', of 47 radio galaxies, and have sum- 
marized our results here. 

(i) The upper limits to the Hi column density towards 
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the cores of large sources presented here range from ~0.5 to 
69xl0 20 cm" 2 . 

(ii) Considering observations of similar sensitivity, we find 
that the detection rate of Hi absorption towards the cores 
of larger (>15 kpc) sources is only ~15 per cent, compared 
with ~57 per cent for CSS and GPS objects, suggesting an 
evolution in the gaseous content of the host galaxies as the 
radio sources age and grow in size. 

(iii) The Hi column density towards the core vs the largest 
linear size plot for the large sources shows no significant 
correlation between these parameters. 

(iv) The distribution of the relative velocity of the princi- 
pal absorbing component towards the cores of large sources 
is not significantly different from that of the CSS and GPS 
objects, although a few of the CSS and GPS objects show 
large blue-shifted velocites (>1000 km s" 1 ) possibly due to 
jet-cloud interactions. 

(v) For the 'cores sample', Hi absorption is detected to- 
wards 4 out of 32 (~13%) FRI objects, and 3 out of 15 
(~20%) FRII sources. Although there have been sugges- 
tions of differences in the torus/disk structure and fuelling 
mechanisms between FR I and FR II sources, with the small 
number of detections there does not appear to be a signifi- 
cant difference in the detection rate between these two types 
of radio sources. 
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Figure Al. Image of J0209+3547 (4C+35.03) with an rms noise 
of 0.2 mjy/beam. The contour levels are 0.6x(-l, 1, 2, 4, 8, 16, 
32, 64, 128) mjy/beam. In all the images the restoring beam is 
shown as an ellipse while the position of the optical host galaxy is 
marked with a (+) sign. The dashed contours represent the nega- 
tive flux density values while solid contours represent the positive 
flux density values. The spectra at the peak intensity pixel (shown 
as (x) sign) for the different components are shown below each 
image and the components have been indicated in the spectra. 
For all the spectra the x- and y-axes represent the heliocentric 
velocity in km s and the normalised intensity respectively. The 
spectra have been normalised using the flux density correspond- 
ing to the peak intensity pixel. 
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Figure A2. Image of J0223+4259 (3C66B) with an rms noise 
of 0.6 mjy/beam. The contour levels are 3x(— 1, 1, 2, 4, 8, 16, 
32, 45) mjy/beam. 
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Figure A3. Image of J0313+4120 with an rms noise of 0.5 
mjy/beam. The contour levels are 2.5x(— 1, 1, 2, 4, 8, 16, 32, 
64, 128) mjy/beam. 
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Figure A4. Image of J0334+3921 (4C+39.12) with an rms noise 
of 0.1 mjy/beam. The contour levels are 0.8x(-l, 1, 2, 4, 8, 16, 
32, 64, 128, 256) mjy/beam. 



APPENDIX A: MAPS AND 21-CM SPECTRA 
OF THE 'CORES SAMPLE' FROM OUR GMRT 
OBSERVATIONS. 
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Figure A5. Image of J0418+3801 (3C 111) with an rms noise of 
0.6 mjy/beam. The contour levels are 6x(-l, 1, 2, 4, 8, 16, 24, 
32, 64, 128) mjy/beam. 
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Figure A7. Image of J0748+5548 (DA 240) with an rms noise 
of 0.37 mjy/beam. The contour levels are 2.40x(-l, 1, 2, 4, 8, 
16, 32, 64, 128, 256) mjy/beam. 
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Figure A6. Image of J0656+4237 (4C 42.22) with an rms noise 
of 0.24 mjy/beam. The contour levels are 0.8x(— 1, 1, 2, 4, 8, 16, 
32, 64, 128) mjy/beam. 
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Figure A8. Image of J0758+3747 (3C 189) with an rms noise 
of 0.5 mjy/beam. The contour levels are 2.0x(— 1, 1, 2, 4, 8, 16, 
32, 64, 128) mjy/beam. 
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Figure A9. Image of J1147+3501 with an rms noise of 0.54 
mjy/beam. The contour levels are 2.057x(-l, 1, 2, 4, 8, 16, 32, 
64, 128, 256) mjy/beam. 
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Figure A10. Image of J1628+3933 with an rms noise of 0.46 
mjy/beam. The contour levels are 5x(-l, 1, 2, 4, 8, 16, 32, 64, 
128) mjy/beam. 
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Figure All. Image of J1744+557 with an rms noise of 0.48 
mjy/beam. The contour levels are 4x(-l, 1, 2, 4, 8, 16, 32, 64, 
128, 256) mjy/beam. 
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Figure A12. Image of J1835+3241 (3C382) with an rms noise 
of 0.5 mjy/beam. The contour levels are 2.35x(-l, 1, 2, 4, 8, 16, 
32, 64, 128, 256) mjy/beam. 
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Figure A13. Image of J1836+1942 (PKS 1834+196) with an 
rms noise of 0.8 mjy/beam. The contour levels are 4.0x(— 1, 1, 
2, 4, 8, 16, 32, 64, 128, 256) mjy/beam. 
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Figure A15. Image of J2338+2701 (3C 465) with an rms noise 
of 0.5 mjy/beam. The contour levels are 2.5x(— 1, 1, 2, 4, 8, 16, 
32, 64, 128) mjy/beam. 




Figure A14. Image of J1842+7946 (3C390.3) with an rms noise 
of 0.59 mjy/beam. The contour levels are 2.35x(— 1, 1, 2, 4, 8, 
16, 32, 64, 128, 256) mjy/beam. 



